Anisotropic Kosterlitz-Thouless Transition Induced by Hard- Wall Boundaries 
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The spatial dependence of the superfiuid density is calculated for the Kosterlitz-Thouless transition in the 
presence of hard-wall boundaries, for the case of a single wall bounding the half-infinite plane, and for a super- 
fluid strip bounded by two walls. The boundaries induce additional vortices that cause the superfiuid density 
to become anisotropic, with the tensor component perpendicular to the wall falling to zero at the wall, whereas 
the component parallel to the wall remains finite. The effects of the boundaries are found to extend over all 
measured length scales, since the correlation length is infinite in the superfiuid phase. 

PACS numbers: 67.40.Hf, 67.40.Kh, 05.10.Cc, 74.78.-w 
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The behavior of a superfiuid near a hard wall is still not well 
understood, even though many experiments are conducted in 
the presence of container walls. The boundary condition at 
the wall is that the superflow normal to the wall should fall 
to zero. There have been many guesses that this can be sat- 
isfied by having the superfiuid density itself fall to zero right 
at the wall, but this has never been verified except perhaps 
in the simplest Landau-Ginsburg models. 1 Dirichlet boundary 
conditions at a wall are often imposed on the order parameter 
in perturbative renormalization expansions, 2 but without any 
real justification, and since the superfiuid density is not sim- 
ply proportional to the order parameter in these theories its 
variation near the wall is often left undetermined. 

Here we calculate directly the superfiuid density near a hard 
wall for the case of the two-dimensional Kosterlitz-Thouless 
(KT) transition, 3 by taking into account the image vortex pairs 
necessary to satisfy the boundary condition. This leads to a 
lowered energy of the pairs in the superfiuid near the wall, 
and a consequent increase in the density of pairs. The results 
show that the superfiuid density does indeed fall to zero at the 
wall, but only for the component measured in the direction 
perpendicular to the wall, while the component parallel to the 
wall remains finite. The presence of the wall thus causes the 
superfiuid density to become an anisotropic tensor, and this 
perturbation is found to extend over all length scales from the 
wall that are being measured. 

We also consider the related problem of two parallel walls 
bounding a long superfiuid strip, a geometry often used in 
experimental measurements. We note that several previous 
studies of the KT transition near hard-wall boundaries have 
been carried out,l but none of these have incorporated the 
anisotropy of the superfiuid density, which we find to be a 
major effect. 

These results should be of interest to experimental stud- 
ies of the KT transition. Although superfiuid 4 He films often 
completely wet their substrates and thus have no boundaries, it 
is now possible to micro-machine channels in silicon to such a 
small thickness that near the bulk lambda transition the helium 
in the channels undergoes a two-dimensional KT transition, 5 
and the geometry and dimensions of the bounding walls can 
easily be varied. The superfiuid density can be measured with 
a small probing flow field in a given direction at frequency ui, 



and the length scale L over which the measurement is made 
is then the vortex diffusion length, 6 which is thought to vary 
as uj^ 1 /' 2 . For helium films at kHz frequencies this length is 
known to be 4-5 orders of magnitude larger than the vortex 
core sizeii 

Consider a vortex pair whose center is a distance z from 
a wall which bounds a semi-infinite half plane of superfiuid 
film. The pair have a separation r and the line joining their 
cores makes an angle 8 with respect to the perpendicular to 
the wall. For convenience we we will scale all of these lengths 
in the problem in units of the vortex core radius ao. Since the 
wall can be replaced by the opposite-sign image vortices, the 
energy of the pair, scaled by ksT, is easily found to be 8 
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where Kq — Ti 2 a° s /m 2 kBT is the bare areal superfiuid den- 
sity in dimensionless form, and E c is the core energy of a 
vortex. Since the energy is a function of the angle with the 
wall, the distribution of the thermally excited vortex pairs 
will be angle-dependent, and this causes the superfiuid den- 
sity to become anisotropic, with tensor components Kn and 
K± . By following the arguments of Machta and Guyes^ for 
the anisotropic KT transition, it can be shown that the scal- 
ing relations for the components of the superfiuid density are 
given by 
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FIG. 1 : Superfluid density components in the half plane a) perpendic- 
ular and b) parallel to the wall, and c) the corresponding fugacities, 
as a function of distance z from the wall. The curves are labeled by 
their values of T/T K t- 



where £ — In r, y is the vortex fugacity, and f(9) equals cos 2 9 
for the parallel moment of the fugacity y\\ and sin 2 9 for y±. 
The lower limit in the integration of Eq. 5 arises when a pair 
near the wall is at the angle 9 m - ln where one of the pair ap- 
proaches the wall to within a core radius, and from the geom- 
etry this minimum value is 9 m i n — arccos((2z — 2)/e l ) when 
2z - 2 < e e , and 6> min = for 2z - 2 > e t . The density 
of pairs n(r, z)with separation between r and r + dr at the 
distance z is given by n(r, z) = (y 2 + y 2 _)/r 2 , and in the limit 
of large z this reduces to the infinite-plane result y 2 /r 2 . 

The recursion relations are iterated using variable-step 
Runge-Kutta techniques, as is the integration of Eq. 5. The it- 



erations start at the bare scale I — and terminate at the scale 
over which the superfluid density is measured, l m = ln(i). 
The initial values at I = are K® = = Kq, and spin- 
wave effects are neglected by assuming a® = a, the liquid 
density. The initial fugacities are y\\ = y± = y = exp(—E c ), 
where we assume the Villain form E c = tt 2 Kq /2. The critical 
value of Kq at Tkt is then Kq c = 0.747583, and the scaled 
temperature is given by T/Tkt = Kq c /Kq. 

Figure 1 shows the components of the superfluid density 
and fugacity as a function of the distance z from the wall for 
the measuring scale i m =Yl (L=1.63xl0 5 ). The superfluid 
density component perpendicular to the wall is zero within 
a core radius of the wall, and then increases smoothly to the 
infinite-plane value for z > L/2. The perpendicular fugacity 
component y± is a maximum at the wall, reflecting the in- 
creased density of pairs parallel to the wall that drives the su- 
perfluid density to zero, and then decreases smoothly with z. 
For distances larger than L/2 it then decreases to the infinite- 
plane value. At high temperatures and very close to the wall 
y± can become larger than one, a high-density regime where 
the recursion relations are not valid, but over most of the 
regime where the superfluid density varies near the wall y± 
is less than one and the results should be accurate. In Figure 
1 both components of the superfluid density vary rapidly near 
Tkt, but there is not a sharp jump to zero because of finite- 
size broadening due to the cutoff of the recursion relations at 
the distance L. This is well known 10 to give a temperature 
broadening above Tkt m the form AT w (6/ InL) 2 , where 
for our parameters the value of the constant b is about 1.7. 

If the distance over which the perpendicular superfluid den- 
sity rises to one-half of its ultimate value at large z is de- 
noted as £1/2, it can be seen in Figure 1 that this quantity 
varies rapidly with temperature. Fits to £]/ 2 over the range 
T/Tkt between 0.75 and 1.0 are best characterized by the 
form£ 1/2 = 0.15£exp(13.3/(l - T/T KT )). This differs 
from the form exp(6/(l — T/Tkt) 1 ^ 2 ) that has been sug- 
gested as a "superfluid" coherence length for the infinite-plane 
casepi but the reason for the difference in exponents is not 
clear; the precise form will require an analytic solution of the 
recursion relations of Eqs. 2-5. 

The component of the superfluid density parallel to the wall 
remains finite at the wall, but with a rapid drop over 10-40 
core radii depending on the temperature. This decrease is due 
to the point in the recursion relations where the pair separa- 
tion is nearly 2z and one of the pairs approaches to within a 
core radius of the wall, effectively canceling with the image 
in the wall. This reduces the energy of the pair considerably, 
resulting in an increased vortex density that lowers the super- 
fluid density. The phase space of such pairs is quite limited, 
however, since larger pairs are cut off at 9 m [ n , and the result- 
ing fugacity component y» is relatively small right at the wall, 
has a tiny peak near 10 core radii for the higher temperatures 
(barely visible in Fig. 1), and then increases slowly with z as 
the phase space increases. For intermediate values of z the 
parallel superfluid density is relatively constant, but then as 
z approaches L/2 it increases sharply as the pair energy is 
reduced only for very small angles 9, and finally is not re- 
duced at all past L/2 where the superfluid density quickly ap- 
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FIG. 2: Superfluid density components a) perpendicular and b) paral- 
lel to the walls of a strip of width W = 1000, and c) the corresponding 
fugacities, as a function of distance z from the left wall. The curves 
are labeled by their values of T /Tkt- 



proaches the infinite-plane value, yu peaks just at L/2 where 
the phase space increases rapidly in conjunction with the low- 
ered pair energy, and then past L/2 it decreases rapidly to the 
infinite-plane value of y. 

If these calculations are repeated at different values of L, 
the results are nearly identical when plotted as a function of 
z/L, with the only differences occurring near Tkt due to the 
finite-!/ broadening of the transition noted above. The effect 
of the boundary extends over all length scales that are being 
measured, a consequence of the infinite Kosterlitz correlation 
length 3 in the superfluid phase below Tkt- The origin of the 
depression of the perpendicular superfluid density at macro- 




FIG. 3: Temperature dependence of the superfluid densities a) per- 
pendicular and b) parallel to the walls of a strip of width W indicated 
by the numbers on the curves, at a distance from the wall z = W/4. 
c) shows the corresponding fugacities, and the dashed portions of the 
curves indicate where the fugacities become larger than one. The 
upper curve in a) and b) is the infinite-plane KT transition. 



scopic values of z can be understood by the excitation of pairs 
oriented parallel to the wall (0 — tt/2), which are primarily 
responsible for the reduction of the perpendicular component. 
It is found from Eq. 1 that the energy to excite such pairs be- 
comes nearly linear in the variable z/r when r is greater than 
z. Hence even at very large z it becomes quite favorable to 
excite pairs with much larger values of r, and these strongly 
reduce the perpendicular superfluid density even at tempera- 
tures well below Tkt- It is only when the maximum value of 
r affecting the superfluid density is limited by the measuring 
length L that the effects of the wall are finally curtailed. 

The case of two hard walls defining a superfluid strip of 
width W has also been investigated using the same techniques, 
with the cutoff at L effectively acting as periodic boundary 
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conditions at the ends of the strip. The image vortices for this 
geometry form two infinite lines of alternating-sign vortices, 
and the bare energy of a pair is given by 

= 2nK [\nr + 
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where k — ir/W, a = r cos 8, and (3 — r sin 8. With this po- 
tential the recursion relations are identical to Eqs. 2-4, while 
in Eq. 5 the the term in the parentheses of the integral is re- 
placed with the argument of the second logarithm in Eq. 6, 
with r = exp(l). Figure 2 shows the calculated superfiuid 
densities and fugacities for a strip of width W = 1000, and 
measuring scale l m — 12. Shown are the values for z be- 
tween and W/2; the results are symmetric about W/2. Near 
z = the characteristics are quite similar to the single-wall 
result, with the perpendicular superfiuid density rising from 
zero at the wall, and the parallel component remaining finite. 
The main difference comes when z approaches W/2. At this 
point both vortices in the pair can approach the walls, greatly 
lowering the energy and causing a sharp dip in both compo- 
nents of the superfiuid density even at very low temperatures. 
An additional effect of the second wall is that the transition 
for the perpendicular component is greatly broadened out, and 
there is also a downward shift in the perpendicular transition 
temperature to well below the infinite-plane Tkt- 

When the calculation is repeated for different values of W, 
the curves are quite similar when plotted versus the scaled 
variable z/W, with the main differences being the broaden- 
ing and T c shift of the perpendicular component. Figure 3 
shows the temperature dependence of the superfiuid densities 
and fugacities at the point z = W/4, for a wide range of W. 
The perpendicular superfiuid density shows the strong down- 
ward shift of T c , which appears to be a logarithmic decrease 



with W, varying as \n(W~ om ). The transition is also greatly 
broadened, with a linear temperature decrease that seems to be 
independent of W. In contrast, the parallel component shows 
about a 15% decrease in T c that is independent of W, and 
then a W^-dependent finite-size temperature broadening vary- 
ing approximately as (b/\nW) 2 , similar to the single-wall 
broadening but with L now replaced by W. 

These main results of anisotropy extending over macro- 
scopic length scales should also apply to superconducting 
films that undergo the KT transition, since strip geometries 
are often used in the measurements. Finite-size effects at the 
KT transition are known to be important in low-T c granular 
films, 11 high-Tc filmsp 2 *^ and Josephson-junction arrays. 14 
However, the superconducting case involves an additional 
length scale, the penetration depth A, 15 and to compare with 
experiments it will be necessary to extend the above theory 
to incorporate that scale. A crude first approximation is to re- 
place the measurement length L in the above calculations with 
A if it is smaller than L. 

It also may be possible to extend the calculation to the 
three-dimensional superfiuid transition 16 by using similar 
techniques on the vortex loops near the wall. This would be 
important for gaining a more complete understanding of the 
critical Casimir effect in helium films, where the boundary 
conditions at the upper and lower film surfaces strongly af- 
fect the magnitude of the Casimir force. 213 The images in the 
wall would again lower the energy to create loops, but unfor- 
tunately analytic solutions for the energy are only known for 
loops parallel to the wall, 17 and numerical techniques will be 
necessary for the energies at other angles. 18 We would expect 
again to find anisotropy in the superfiuid density, but since in 
3D the correlation length is finite this would only occur within 
a correlation length from the wall. 
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